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ABSTRACT 

An ever-increasing body of evidence suggests that weakly interacting massive parti- 
cles (WIMPs) constitute the bulk of the matter in the Universe. Experimental data, 
dimensional analysis and Standard Model particle physics are sufficient to evaluate 
and compare the performance of detectors searching for such particles either directly 
(e.g. by their scattering in germanium detectors), or indirectly (e.g. by observing 
their annihilation into neutrinos in underground detectors). We conclude that the 
direct method is superior if the WIMP interacts coherently and its mass is lower 
or comparable to the weak boson mass. In all other cases, i.e. for relatively heavy 
WIMPs and for WIMPs interacting incoherently, the indirect method will be com- 
petitive or superior, but it is, of course, held hostage to the successful deployment 
of high energy neutrino telescopes with effective area in the ~ 10 4 -10 5 m 2 range 
and with appropriately low threshold. The rule of thumb is that a kilogram of ger- 
manium is roughly equivalent to a 10 4 m 2 neutrino telescope, although the signal- 
to-noise is, at least theoretically, superior for the neutrino detector. The energy 
resolution of the neutrino telescope may be exploited to measure the WIMP mass 
and suppress the background. A kilometer-size detector probes WIMP masses up 
to the TeV-range, beyond which they are excluded by cosmological considerations. 



1. Introduction and Results 



It has become widely accepted that most of our Universe is made of cold dark 
matter particles. Big bang cosmology implies that these particles have interactions 
of order the weak scale, i.e. they are WIMPs! We briefly review the argument which 
is sketched in Fig. 1. In the early Universe WIMPs are in equilibrium with photons. 
When the Universe cools to temperatures well below the mass m x of the WIMP their 
density is Boltzmann-suppressed as exp(— m x /T) and would, today, be exponentially 
small if it were not for the expansion of the Universe. At some point, as a result of this 
expansion, WIMPs drop out of equilibrium with other particles and a relic abundance 
persists. The mechanism is analogous to nucleosynthesis where the density of helium 
and other elements is determined by competition between the rate of nuclear reactions 
and the expansion of the Universe. 



At high temperatures WIMPS are abundant and they rapidly convert into lighter 
particles. Also, as long as they are in equilibrium, lighter particles interact and 
create WIMPs. The situation changes rapidly after the temperature drops below the 
threshold for creating WIMPs, T < m x . The WIMP density falls exponentially as a 
result of their annihilation into lighter particles. When the expansion of the Universe 
has reduced their density to the point where annihilation is no longer possible, a relic 
density "freezes out" which determines the abundance of WIMPs today. This density 
is just determined by the annihilation cross section; for a larger cross section freeze- 
out is delayed resulting in a lower abundance today and vice versa. The scenario is 
sketched in Fig. 1 where the density of WIMPs (in the comoving frame) is shown as 
a function of time parametrized as the inverse of the temperature m x /T. 
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Fig. 1 

For WIMPs to make up a large fraction of the Universe today, i.e. a large fraction 
of Q, their annihilation cross section has to be "just right". The annihilation cross 
section can be dimensionally written as a 2 /m x , where a is the fine-structure constant. 
It then follows that 

Q oc 1/cr oc m x . (1) 

The critical point is that for Q ~ 1 we find that m x ~ m^, the mass of the weak 
intermediate boson. There is a deep connection between critical cosmological density 
and the weak scale. Weakly interacting particles which constitute the bulk of the mass 
of the Universe remain to be discovered. It may not be an accident that the unruly 
behavior of radiative corrections in the Standard Model also requires the existence of 
such (supersymmetric?) particles. 
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When our galaxy was formed the cold dark matter inevitably clustered with the 
luminous matter to form a sizeable fraction of the 



p x = 0.4 GeV/cm 3 



(2) 



galactic matter density implied by observed rotation curves. Unlike the baryons, the 
dissipationless WIMPs fill the galactic halo which is believed to be an isothermal 
sphere of WIMPs with average velocity 



In summary, we know everything about these particles (except whether they really 
exist!). We know that their mass is of order of the weak boson mass; we know that 
they interact weakly. We also know their density and average velocity in our Galaxy 
given the assumption that they constitute the dominant component of the density of 
our galactic halo as measured by rotation curves. 

For a first look at the experimental problem of how to detect these particles it is 
sufficient to recall that they are weakly interacting with masses in the range 



WIMPs have a mass of order the weak boson mass, in the tens of GeV to several TeV 
range. Lower masses are excluded by accelerator and (in) direct searches with existing 
detectors while masses beyond several TeV are excluded by cosmological considera- 
tions. Two general techniques, referred to as direct (D) and indirect (ID), are pursued 
to demonstrate the existence of WIMPs.! In direct detectors one observes the energy 
deposited when WIMPs elastically scatter off nuclei. The indirect method infers the 
existence of WIMPs from observation of their annihilation products. WIMPs will 
annihilate into neutrinos; massive WIMPs will annihilate into high-energy neutrinos 
which can be detected in high-energy neutrino telescopes.! Throughout this paper we 
will assume that such neutrinos are detected in a generic Cherenkov detector which 
measures the direction and, to some extent, the energy of a secondary muon pro- 
duced by a neutrino of WIMP origin in or near the instrument. It can oalso detect 
the showers initiated by electron-neutrinos. 

The indirect detection is greatly facilitated by the fact that the sun represents 
a dense and nearby source of accumulated cold dark matter particles.! Galactic 
WIMPs, scattering off nuclei in the sun, lose energy. They may fall below escape 
velocity and be gravitationally trapped. Trapped WIMPs eventually come to equi- 
librium temperature and accumulate near the center of the sun. While the WIMP 
density builds up, their annihilation rate into lighter particles increases until equi- 
librium is achieved where the annihilation rate equals half of the capture rate. The 
sun has thus become a reservoir of WIMPs which we expect to annihilate mostly into 
heavy quarks and, for the heavier WIMPs, into weak bosons. The leptonic decays of 
the heavy quark and weak boson annihilation products turn the sun into a source of 



v x = 300 km/sec . 



(3) 



tens of GeV < m v < several TeV . 



(4) 



3 



high-energy neutrinos with energies in the GeV to TeV range, rather than in the keV 
to MeV range typical for neutrinos from thermonuclear burning. 

The performance of future detectors is determined by the rate of elastic scattering 
of WIMPs in a low-background, germanium detector and, for the indirect method, 
by the flux of solar neutrinos of WIMP origin. Both are a function of WIMP mass 
and of their elastic cross section on nucleons. In standard cosmology WIMP capture 
and annihilation interactions are weak, and we will suggest that, given this constraint, 
dimensional analysis is sufficient to compute the scattering rates in germanium detec- 
tors as well as the neutrino flux from the measured WIMP density in our galactic halo. 
We derive and compare rates for direct and indirect detection of weakly interacting 
particles with mass m x ~ assuming 

1. that WIMPs represent the major fraction of the measured halo density, i.e. 

0.4 GeV in6 cm 1.2 x 10 7 _ 2 _ x 
y = n y v y = 3 x 1(T — = cm z s , (5) 

m x CUT S T^xGeV 

where m X GeV = GeV) is in GeV units. 

2. a WIMP-nucleon interaction cross section based on dimensional analysis 

a{ X N) = (G F m 2 N ) 2 — = a DA = 6 x 1(T 42 chi 2 . (6) 
m w 

3. that WIMPs annihilate 10% of the time in neutrinos (this is just the leptonic 
branching ratio of the final state particles in the dominant annihilation channels 
XX ~^ W + W~ or QQ, where Q is a heavy quark). 

Clearly the cross section for the interaction of WIMPs with matter is uncertain. 
Arguments can be invoked to raise or decrease it. Important points are that i) our 
choice represents a typical intermediate value, ii) all our results for event rates scale 
linearly in the cross section and can be easily reinterpreted, and iii) the comparison 
of direct and indirect event rates is independent of the choice. 

Our conclusions will not be surprising.^ We find that the direct method is superior 
if the WIMP interacts coherently and, if its mass is lower or comparable to the weak 
boson mass mw- In all other cases, i.e. for relatively heavy WIMPs and for all 
WIMPs interacting incoherently, the indirect method is competitive or superior, but 
it is, of course, held hostage to the successful deployment of high energy neutrino 
telescopes with effective area in the ~ 10 4 -10 6 m 2 range and with appropriately low 
threshold. Especially for heavier WIMPs the indirect technique is powerful because 
underground high energy neutrino detectors have been optimized to be sensitive in 
the energy region where the neutrino interaction cross section and the range of the 
muon are large. A kilometer-size detector probes WIMP masses up to the TeV-range, 
beyond which they are excluded by cosmological considerations. 



4 



For high energy neutrinos the muon and neutrino are aligned, with good angular 
resolution, along a direction pointing back to the sun. The number of background 
events of atmospheric neutrino origin in the pixel containing the signal will be small. 
The angular spread of secondary muons from neutrinos coming from the direction of 



the sun is well described by the relatione ~ 1.2y ^ E^TeV). Measurement of muon 
energy, which may be only up to order of magnitude accuracy in some experiments, 
can be used to infer the WIMP mass from the angular spread of the signal. The 
spread contains information on the neutrino energy and, therefore, the WIMP mass. 
More realistically, measurement of the muon energy can be used to reduce the search 
window around the sun, resulting in a reduced background. 

Our analysis will quantify all statements above in a simple and totally transparent 
framework. It finesses all detailed dynamics and gives answers that are sufficiently 
accurate considering that the mass of the particle has not been pinned down. 

Before proceeding, we comment on our ansatz for the elastic WIMP-nucleon scat- 
tering cross section. The simplest dimensional analysis implies that the cross section 
is G 2 F m 2 N . This correctly describes the Z-exchange diagram of Fig 2a, which is of the 
form 

r 2 m2 N m l m 

(m N + m x y 

For coherent interactions, which we will emphasize throughout this paper, there 
is an additional suppression factor associated with the exchange of the Higgs particle 
with a mass of order of the weak boson mass; see Fig 2b. In the specific diagram 
shown the Higgs interacts with the heavy quarks in the gluon condensate associated 
with the nucleon target. It is of the form 

„ ~ G F9 % « 



(mjv + rnx) 2 m 



where gu ~ VGf ttin describes the condensate. Conservatively, we will use the sup- 
pressed WIMP interaction cross section which is appropriate for coherent scattering. 



(a) 

N 



(b) 1 H 



Q 

Fig. 2. Examples of (a) incoherent and (b) coherent WIMP-nucleon interactions. In (b) the gluon 
is a constituent of the target nucleon and Q is a heavy quark. 



2. Indirect Detection (ID) 

The number of solar neutrinos of WIMP origin can be calculated in 5 easy steps 
by determining 
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• the capture cross section in the sun, which is given by the product of the number 
of target nucleons in the sun and the elastic scattering cross section 

a Q = f [1.2 x 10 57 ] a DA . (9) 

This includes a focussing factor / given, as usual, by the ratio of kinetic and 
potential energy of the WIMP near the sun. It enhances the capture rate by a 
factor 10. 

• the WIMP flux from the sun which is given by 

= <P x a & /And 2 , (10) 
where d = 1 a.u. = 1.5 x 10 13 cm. 

• the actual neutrino flux, which is obtained after inclusion of the branching ratio. 
From (§),(§ and ®,(0) 

3 x 10~ 5 

(j) u = KT 1 x (j) Q = cm" 2 s- 1 . (11) 

I^X GeV 

• the probability to detect the neutrinoj which is proportional to 

P = pa u R fl , with 

p = Avagadro # = 6x 10 23 

cr„ = neutrino interaction cross section = 0.5 x 10~ 38 E u (GeV) cm 2 
Rfj, = muon range = 500 cm £7^ (GeV) 
or 

P = 2xlO- 13 m 2 GeV (12) 
Here we assumed the kinematics of the decay chain 

xx - w + w- 

with E v = \m x (this would be \m x for Q decay) and E^ = \E V = 

• finally, dN 1D /dA = <p u P = 1.8 x 10~ 6 m xGeV (year)" 1 (m 2 )" 1 (13) 
where dNis/dA represents the number of events from the sun per unit area (m 2 ) 
detected by a neutrino telescope. 

The linear rise of a u , with energy, which are the origin of the good detection 
capability of neutrino telescopes for large WIMP masses, are valid approximations 
up to 



m x m 2 w 
E u ^ -7T £ ' and 

HI 

E„ ^ -f > 500 GeV, 
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so the approximations are valid for m x well into the TeV mass range. This is sufficient 
as m x ^> 1 TeV is cosmologically unacceptable. 

3. Direct Detection 

The event rate in a direct detector is proportional to the WIMP cross section, 
flux and the density of targets m^ 1 , i.e. 

dm rriN 

where represents the number of direct events per unit of target mass. 

We can now summarize our results so far by comparing a 10 4 m 2 neutrino detector, 
an area typical of the instruments now being deployed, with a kilogram of hydrogen: 

dN m /dA = 1.8 x l(r 2 m X GeV (10 4 m^-^year)" 1 
dN D /dM = L4 (kg)" 1 (year)" 1 

^xGcV 

dN D /dM / lQ 4 m 2 \ 7.8 x 10 1 

dN m /dA [~^r) ~ m 2 GeV ( } 

Direct detection is superior only in the mass range m x < 10 GeV, but this region is, 
arguably, ruled out by previous searches. Indirect detection is the preferred technique. 
This straightforward conclusion may, however, be invalidated when WIMPs interact 
coherently and targets other than hydrogen are considered. We discuss this next. 

4. Coherent Nuclear Enhancements 

The nuclear dependence of the event rates resides in 

• the target density factor in Eq. (14). The mass of the target nucleus 
rriA = Am N is substituted for m N . 

• the coherent enhancement factor "A 2 ", 

• the nuclear dependence of the cross section is obtained by the substitution 
incoherent 

, 2 m 2 N m 2 x {Am N fm 2 x 



F {rriN + m x ) 2 F (ArriN + m x ) 2 
coherent 

~ „ „ GWh m2 » m i * ^ Gf (i^y , ^ m "i 2 %, 

(m N + m x y ?t% \m w J (Am N + m x ) 2 
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The coherent enhancement factor for a nucleus A, including a factor A 1 for 
the target density, is given by 

1 A 2 (Am N ) 2 m 2 x (m N + m x ) 2 _ ^ (m N + m x ) 2 
A{Am^ + m x ) 2 m 2 v m x (Am^ + m x ) 2 

Below we list a sample of enhancement factors having in mind germanium instead 
of hydrogen for direct detection and oxygen or iron nuclei for capture in the sun. 
Their solar abundance relative to hydrogen are 1.1 and 0.2 percent, respectively. 



Table 1. Nuclear Enhancement Factors 





m x = 50 


m x = 500 


m x = 2000 


Ge (A = 


76) 


7.7 x 10 4 


3.9 x 10 5 


4.2 x 10 5 


(A = 


16) 


2.5 x 10 3 


3.9 x 10 3 


4.1 x 10 3 


Fe (A = 


56) 


4.3 x 10 4 


1.4 x 10 5 


1.7 x 10 5 



Because of the complications associated with nuclear form factors, calculations 
considering just oxygen or iron capture in the sun bracket detailed computations.0 

5. Event Rates for WIMPs with Coherent Interactions 

After inclusion of above coherence factors in Eq. ( |T5"D we obtain the data rates 
listed below: 



Table 2. Event Rates per year for Direct and Indirect Detection 





ID/10 4 m 2 D/kg D/ID ratio (kg/10 4 m 2 ) 


m x = 50 


(H) 9 x 10" 1 (H) 2.8 x 10~ 2 3.1 x 10" 2 
(O) 2.5 x 10 1 (Ge) 2.2 x 10 3 8.8 x 10 1 


m x = 500 


(H) 9 (H) 2.8 x 10~ 3 3 x 10" 4 
(O) 4 x TO 2 (Ge) 1.1 x 10 3 2.7 


m x = 2000 


(H) 3.6 x 10 1 (H) 7 x 10~ 4 2 x 10~ 5 
(O) 1.6 x 10 3 (Ge) 2.9 x 10 2 1.8 x 10" 1 



We here assumed a direct detector made of germanium. Conservatively, only 
oxygen, the dominant element in the solar capture rate, was considered in calculating 
the indirect rates. We chose 50 GeV and 2 TeV WIMP as illustrative masses in order 
to bracket the appropriate range with an illustrative central value of 500 GeV. One 
should keep in mind that a 500 GeV WIMP is well out of reach of present as well as 
future accelerator searches. 



1 _i_ UhL 

m N 



m N 



A + 



(16) 



8 



The ratio of direct to indirect events, which is independent of the WIMP-nucleon 
cross section, can be summarised by the following equation: 



D _ dN D /dM 
ID ~ dN 1D /dA 



7.8 x 10 1 



N(A B ) 



m 



xGeV 



N(A m ) \p{A m )/p{H)\ 



(17) 



with 



N(A) = A 3 



I _|_ m x 

m N 



As in Eq. (|15D the units are 



10 a m z 
kg 



Ad and Aid are the atomic numbers appropriate 



for the nuclei involved in the direct detection and capture in the sun, respectively. 
The latter is weighted by its relative mass abundance \p(Am)l p{H)] in the sun and 
a summation over elements is understood. This relative evaluation of the two ex- 
perimental techniques is in very good qualitative agreement with a similar analysis 
performed in the context of supersymmetry.i 



6. Final Event Rates 

Our simple evaluations, made so far, overestimate the indirect rates for very heavy 
WIMPS because high energy neutrinos, created by annihilation near the core, may 
be absorbed in the sun. Absorption is stronger for neutrinos and, therfore, mostly 
antineutrinos form the signature for very heavy WIMPS. The probability that an 
antineutrino escapes without absorption is well parametrized by (1 + 3.8 x 10~ 4 E U )~ 7 , 
where E v ~ m x /2. The final rates for indirect detection are 



dN lD /dA ~ {l.8 x 10 2 m xGc y\ < 0.01L4 2 





"l + 


m x 


2' 


jo.oiiA 3 


mjv 


> 




A + 


m x 






m N _ 





>{l + 1.9x l(T 4 m X Gev} 7 



(19) 

The relative merits of the two methods are summarised in the following table, 
which establishes that a kilogram of germanium and a 10 4 m 2 are competitive. 



Table 3. Event rates and signal to background (Ay B). 



m x (GeV) 


Direct (/kg/year) 


Indirect (/10 4 


m 2 /year) 




events N/ B 


events 


N/B 


50 


2.2 x 10 3 7 


2.3 x 10 1 


~ 1 


500 


1.1 x 10 3 7 


2 x 10 2 


~ 10 2 


2000 


2.9 x 10 2 1 


1.7 x 10 2 


~ 10 4 



At the lower energy the event rates for the indirect method are underestimated 
because also the Earth is a source of neutrinos of WIMP origin. 
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We conclude that the direct method yields more events for the lower masses, even 
when compared to a 10 6 m 2 detector. As expected, the indirect method is competitive 
for heavier WIMPs with a detection rate growing like or m 2 . A 10 5 m 2 covers the 
full WIMP mass range, even if the WIMPs do not coherently interact with nuclei in 
the sun; see Table 2. These conclusions are reinforced after considering the signal- 
to-noise for both techniques. We show this in the next section. A summary of our 
results is shown in Fig. 3. 



#events 
/kg/year 



Direct 




10 3 5000 



#events 
/10 4 m 2 /year 



Indirect 




10 3 5000 



111-, 



nu 




Direct: 1 kg Ge 



Indirect: 10 4 m 2 telescope 

absorption of v in sun 



o(xiV) = 6xlO 42 cm 2 



Fig. 3. Summary of the results. The results shown are for a(xN) = 6 x 10 42 cm 2 . All event rates 
scale linearly in a{xN). The relative direct and indirect rates are independent of a{xN). 
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7. Backgrounds 



Indirect Background. For the indirect detection the background event rate is deter- 
mined by the flux of atmospheric neutrinos in the detector coming from a pixel around 
the sun.B The number of events in a 10 4 m 2 detector is ~ 10 2 /£' /i (TeV) and the pixel 

size is determined by the angle between muon and neutrino ~ 1.2°/ ^jE^TeV). Using 
the kinematics ~ 



B 



ID 



m x /4 we obtain 
10 2 /^(TeV) 



1.1 x 10 5 



2tt, 



1.2° 



per 10 4 m 2 per year 



This is only valid for large m x , 
imation we obtain 



i.e. 



for E^ = 



m X GeV 

m x /4 > 100 GeV. Without this approx- 



Table 4. 



E^GeV) 


# b k S,d- 4 ev § nts 
m 10 m 

in 2n 


# pixels of solar 
size in 2n 


bkgd events 
per 10 m 
per pixel, per year 


10 


3200 


140 


23 


100 


1060 


1.4 x 10 3 


0.8 


1000 


110 


1.4 x 10 4 


8 x lO" 3 



For large m x the signal to background ratio is 

N\ dNi D /dA q 

Clearly, the extremely optimistic predictions for signal-to-noise are unlikely to survive 
the realities of experimental physics. One expects, typically, to measure muon energy 
only to order-of-magnitude accuracy in the initial experiments. The energy of showers 
initiated by electron neutrinos should be determined to a factor 2. It is not excluded 
that future, dedicated experiments may do better. The conclusion that high energy 
muons pointing at the sun represents a superb signature, is unlikely to be invalidated. 

Direct Background: about 300 events per year per kg.! Signal-to-noise therefore 
exceeds unity up to 2 TeV WIMP mass. 

These considerations were used to estimate the signal-to-noise N/B in Table 2. 

8. Dynamics? 

We emphasize that above considerations are valid for the specific and much stud- 
ied example where the lightest supersymmetric particle is Nature's WIMP.Q Clearly 
dynamics, which is now defined, can alter our conclusions, but only in "conspirato- 
rial" ways. In the favored scenario the WIMP is the stable neutralino, i.e. the lightest 
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state composed of the supersymmetric partners of the photon, neutral weak boson 
and the two Higgs particles: 



x = z u W 3 + z 12 B + zwHx + z u H 2 ■ 



(20) 



In this specific case of supersymmetry the diagrams in Fig. 2a,b are proportional to 



where qa is the axial coupling of the nucleon and the weak angle. Even though (22) 
is not completely general, it is clear that the rates can be suppressed in a scheme where 
the neutralino is mostly gaugino-like, i.e. z± 3 = z u = 0. Although such scenarios have 
been suggested, there is no consensus. Dynamics can, on the other hand, increase 
rates as well, sometimes by well over an order of magnitude, over and above the 
rates obtained from dimensional analysis in this paper. Our qualitative conclusions 
are valid, at least in some average sense, in supersymmetry. Our results do, in fact, 
closely trace the supersymmetry prediction of reference 1 for the choice of Higgs 
coupling an = 1, in their notation. 

We feel that the development of detectors should be guided by an analysis like 
ours rather than by dynamics of theories beyond the standard model for which there 
is, at present, no experimental guidance. 
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